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METHOD OF PREDICTING FORMATION TEMPERATURE 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates generally to the prediction of formation temperatures in a 
subsurface formation and, more particularly, to the prediction of the temperature of a 
hydrocarbon formation temperature. 

2. Description of Related Art 

Hydrocarbon fluids, such as oil and natural gas, are obtained from a subterranean 
geologic formation, referred to as a reservoir, by drilling a well that penetrates the hydrocarbon- 
bearing formation. An understanding of the undisturbed reservoir temperature is desirable for 
numerous applications involved in the drilling, completion and production of reservoir fluids. 
These applications may include, for example: drilling fluid and cement slurry design: log 
interpretations; corrosion tendencies in wellbore tubulars and downhole equipment; hydrocarbon 
reserve estimation; flow assurance design; and estimations of geothermal energy, etc. 

hi drilling operations, the formation temperature has a direct bearing on drilling fluid 
rheology and therefore has to be considered in drilling fluid and wellbore design. The formation 
temperature directly impacts cement dehydration and cure times, and therefore needs to be 
considered in the design of casing and cementing programs. The interpretation of electric logs 
requires accurate formation resistivities, which are dependent on temperature. In production and 
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well-control operations, accurate computations of fluid flow rates are important. Fluid 
temperature, both as a function of depth and elapsed time, dictates fluid properties such as 
density and \ iscosity, and therefore influences the pressure drops and/or the ma.xinium allowable 
production rates that can be achieved. Flow assurance design considerations, such as hydrate 
5 formation and paraffin deposition prevention, depend on an accurate knowledge of the reservoir 
temperature. 

As a wellbore is drilled, a temperature disturbance is introduced by the circulating drilling 
fluids, thereby cooling the formation around the borehole. The initial undisturbed formation 
temperature exists only at a certain distance away from the wellbore. During the circulation of 

10 fluids, often referred to as "mud", the temperature within the borehole drops and reaches a 
pseudo-steady state condition in a very short time. After a certam period of time, the temperature 
within the wellbore during fluid circulation can be considered constant. Earlier studies have 
indicated that a constant temperature difference between the bottom-hole fluid and wellbore wall 
is achieved almost immediately and maintained throughout the life of the wellbore fluid 

15 circulation process. This means that the heat transmission from the formation to the wellbore 
during wellbore fluid circulation is a constant heat flux dominated process. Therefore, during 
mud circulation, the heat transfer between the wellbore and the formation can be described with a 
constant heat flux solution of an infinite reservoir. See Raymond. L.R.: "Temperature 
Distribution in a Circulating Drilling TluiJ". JPT. March 1969; ami Schoeppcl. R.J.. Bennet, 

20 R.E.: "Numerical Simulation of Borehole and Formation Temperature Distrihutions While 
Drilling to Total Depth ". SPE paper 3364. presented October 1971. 
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The amount of departure from the undisturbed formation temperature during drilling and 
completion operations depends upon several factors, such as, the original temperature 
distribution, the physical properties of the reservoir rock and the drilling completion fluids. 
Fluid circulation rates and duration, and the tubular and cementation design used on the well are 
5 also factors that influence the temperature profile. Formation temperatures are often estunated 
by using temperature measurements taken inside the wellbore, often in conjunction with well 
logging and fluid sampling. 

The process of formation fiuid .sampling typically involves the lowering of a sampling 
tool into the wellbore. The sampling tool collects one or more samples of formation fluid by the 

JO engagement between a probe module of the sampling tool and the wall of the wellbore. 
Embodiments of samphng tools may comprise more than a single probe, such as with dual-probe 
or multi-probe modules, enabling the sampling of differing sites within the formation within a 
single deployment of the sampling tool. There are several commercially available sample tools 
available, for example the Modular Dynamics Formation Tester (MDT'^') made by 

15 Schlumberger. the Reservoir Characterization Instrument (RCT^') from Baker Atlas, and the 
Reservoir Description Tool (RDT'^') tool made by Halliburton. 

The Modular Dynamics Formation Tester (MDT) formation testing tool, owned and 
provided by Schlumberger operates by creating a pressure differential across an engagement of a 
probe module with the wellbore to induce formation fluid flow mto one or more sample 
20 chambers within the sampling tool. 1 his and similar processes are described in U.S. Patents Nos. 
4,860,581: 4.936.139 (both assigned to Schlumberger). Due to the changes in the temperature 
field surrounding the wellbore discussed above, the temperature data acquired by the MDT is 
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typically lower than the actual static formation temperature, because of short sampling time. One 
distinct feature of wireline pretest and sampling is that the flow regime is primarily controlled by 
three-dimensional (3-D) spherical or radial tlow where the probe functions as a point sink. 
Therefore, the specific difficulty in determining the original formation temperature is the 
5 calculation of the fluid temperature at the probe during the recording, which is associated with 3- 
D spherical flow. 

Another sampling tool is the Reservoir Characterization Instrument (RCI), provided by 
Baker Atlas. It can comprise an optical analyzer, named SampleView"^^', that can be used to 
monitor contamination levels within sample formation fluids pumped through the tool, and can 

10 be run with other reservoir characterization sensors. Examples of other reservoir characterization 
sensors include pressure sensors and sensors that measure the apparent dielectric constant of the 
sample fluid, thereby distinguishing oil, gas and water within the sample fluids. The quality of a 
reservoir fluid sample and the time required to acquire the sample can be predicted utilizing a 
three-dimensional fluid flow simulation model and input data acquired from the RCI, such as 

|5 formation pressure, formation permeabilities, and formation fluid properties. 

Still another sampling tool that can be utilized with the present invention includes the 
Reservoir Description Tool (RDT) tool manufactured by Halliburton. It can comprise a modular 
apparatus that uses nuclear magnetic resonance (NMR) techniques for making downhole NMR 
measurements of the formation fluid samples, as described in U.S. Patent 6,1 1 1.408 to Blades ct 
20 al. 
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Referring to the attached drawings. Figure 1 ilkistrates a prior art representative 
drilling/production platform 10 having a tubular string 12 extending into a wellbore 14. The 
wellbore 14 has penetrated subterranean formations 16, and intersects a productive reservoir 18. 
A casing string 20 lines the well and provides support and isolation of the wellbore 14 from the 
formations 16 and bodies of water 22. Wellbore drilling or completion fluids 24, commonly 
referred to as "mud", are typically circulated down the tubular string 12 and up the wellbore 14. 
The circulation of wellbore fluids 24 results in the cooling of the reservoir 18 around the 
wellbore 14. Upon the cessation of fluid circulation, the tubing 12 can be removed from the 
wellbore 14. 

Figure 2 illustrates a prior art representative drilling/production platform 10 having a 
downhole tool 26 inserted into the wellbore 14 on a wireline 28. The downhole tool 26 can 
comprise a formation-testing tool capable of collecting one or more samples of formation fluid, 
such as, for example, the Modular Dynamics Formation Tester (MDT) formation-testing tool. In 
addition to obtaining a sample of the formation fluid coming from the reservoir 1 8. the downhole 
tool 26 may also collect data such as temperature and pressure readings. Embodiments of the 
downhole tool 26 can be run into the well on a tubing string, slickline, wireline, or by other 
means of positioning the tool within the reservoir 18. The temperature T of the formation fluids 
from the reservoir adjacent the wellbore is lower than the original undisturbed temperature of the 
formation Te, sometimes referred to as the initial or static reservoir temperature or formation 
temperature. When the circulation of the wellbore fluids 24 is stopped, the wellbore fluid 24 
temperature Tm begins to increase due to the influence of the higher temperatures within the 
reservoir 1 8. 
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Figure 3 illustrates a prior art embodiment of an MDT formation testing tool 30 that can 
be utilized in both formation fluid sampling and pretest operations. Various pretest operations 
can invok e flowing formation fluids for a desired amount of time to obtain a specific quantity of 
fluid removal from the formation, or can comprise flowing a well until a desired pressure 
5 drawdown has occurred, in order to measure the pressure recovery or buildup rate. The tool 30 
comprises a pump 32 that induces flow of reservoir fluid into a fluid probe 46, along a flowiine 
34, and then out to the wellbore 14. The flow of reservoir fluid can last an extended period of 
time to reduce contamination and obtain a better quality reservoir fluid sample, or to provide a 
desired amount of pretest flow prior to a pressure buildup test or other formation analysis. If it is 

10 desired to collect a sample, the opening of seal valves 36, 38 divert a portion of the fluid flow 
into a sample chamber 40. Any initial fluid that is in the sample chamber 40 can be flushed with 
the formation fluid. A piston 42 within the tool 30 can move and displace fluid from a buffer 
chamber 44, thus allowing formation fluid to enter the sample chamber 40. The closing of the 
seal valves 36. 38 contain the fluid sample within the sample chamber 40 for removal from the 

15 wellbore 14 and analysis. 

The ability to flow the formation fluid prior to taking a sample enables the MDT tool 30 
to provide a more representative sample of the formation fluid by minimizing near wellbore 
factors such as lost drilling fluids and residual drilling mud from contaminating the sample 
fluids. The MDT tool 30 can comprise a temperature sensor that records the temperature of the 
20 reservoir fluids passing through the flowiine 34 during the testing time period. 

Current methods to determine the initial reservoir temperature are typically based on 
extrapolated shut-in temperature recordings. These methods typically require long shut-in 
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periods and result in estimates that are lower than the true reservoir temperature. Complete 
temperature recovery in the area near the wellbore may take anywhere from a few hours to a few 
months, depending on the formation, well characteristics, and the mud circulating time. Since a 
long waiting period for complete temperature recovery can result in a significant increa.se in 
drilling costs; a less time consuming method is needed to calculate static reservoir temperature 
using early shut-in and test data. 



SUMMARY OF THE INVENTION 

10 One embodiment of the present invention is a method of calculating a static formation 

temperature in a reservoir penetrated by a wellbore. The method comprises estimating the static 
formation temperature and calculating a formation fluid temperature at the wellbore. the 
calculation based, in part, on the estimated static formation temperature. The temperature of a 
sample of formation fluid at the wellbore is measured. The calculated formation fluid 

15 temperature at the wellbore is compared with the measured temperature of the sample of 
formation fluid. The static formation temperature is predicted by altering the estimate of the 
static formation temperature until an error between the calculated formation fluid temperature at 
the wellbore and the mea.sured fluid formation temperature is minimized. 

20 An alternate embodiment of the present invention is a method of calculating a static 

formation temperature in a reservoir penetrated by a wellbore comprising: estimating the static 
formation temperature in the reservoir and a wellbore fluid temperature. A calculated formation 
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fluid temperature at the wellbore versus time profile is created for fluid removed from tlie 
formation by a sink probe, based upon, in part, the estimates of the static formation temperature 
in the reserv oir and the vvcUbore fluid temperature. The temperature of the formation fluid at the 
wellbore removed from the formation by the sink probe is measured, and a measured fluid 
5 formation temperature at the wellbore versus time profile is created. The measured fluid 
formation temperature at the wellbore versus tmie profile is compared to the calculated formation 
fluid temperature at the wellbore versus time profile. The static formation temperature is 
predicted by altering the estimates of the static formation temperature in the reservoir and a 
wellbore fluid temperature until the enor between the measured fluid formation temperature at 
10 the wellbore versus time profile to the calculated fonnation fluid temperature at the wellbore 
versus time profile is minimized. 

BRIEF DESC RIPTION OF THE DRAWINGS 

The invention may be understood by reference to the following description taken in 
conjunction with the accompanying drawings, in which like reference numerals identify like 
elements, and in which: 

Figure 1 illustrates a prior art representative drilling production platforni; 

Figure 2 illustrates a prior art representative drilling/production platform comprising a 
20 formation-testing tool; 
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Figure 3 schematically illustrates a prior art embodiment of an MD I formation testing 

took 

Figure 4 graphically illustrates an example of radial temperature distribution after 
wellbore fluid circulation has ceased: 

5 Figure 5 graphically illustrates an example of calculated three-dmiensional fluid flow 

profdes at various drawdown times. A spherical model is drawn within Cartesian coordinates for 
illustration purposes; 

Figure 6 graphically illustrates an example of calculated and measured temperatures at a 
sink probe versus elapsed time after wellbore fluid circulation has ceased; and 

10 Figure 7 is a schematic process diagram of the methodology of one embodmient of the 

present invention. 

DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

15 Temperature distribution prediction 

During the circulation of fluids within a wellbore, the temperature within the wellbore 
drops and reaches a pseudo-steady .state condition in a very short time. The temperature within 
the wellbore during fluid circulation can be considered constant. A constant temperature 
difference between the bottom-hole tluid and wellbore w all is achieved almost immediately and 
20 maintained throughout the life of the wellbore fluid circulation process. This means that the heat 
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transmission between the wellbore and formation is a constant-heat-flux dommated process. 
Therefore, during mud circulation, the heat transfer between the wellbore and formation can be 

described v\ ith a constant-heat-tlux solution of an intlnite reservoir. 

Applying the constant-heat-flux hnear source solution, the temperature f at a point m the 
5 formation at any time t can be expressed as: 

r^Te--^- t(rj) (1) 



f(rj) = --E,(-'^) (2) 
2 4ar 

Where Te is original formation temperature, Q is heat flow rate (= 27rrHq), k| is thermal 
10 conductivity of formation, H is the reservoir thickness, q is heat flux, E, is the exponential 
integral function, a is the thermal diffusivity of the formation, r is the radius from the wellbore. 
and t is the tune for mud circulation. 

Based on the energy balance of the overall system, the heat rate can be written as: 



27tH{Te-rm) 



Where hr is the heat transfer coefficient for fluid within the wellbore and r.^ is the 
wellbore radius. Tm is the wellbore fluid temperature within the wellbore at the time of 
circulation cessation, which is a time-dependent variable. Since Tm is usually not available in 
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practice, it is treated as an adjustable parameter. Within this apphcation the term wellbore fluid 
refers to the fluids that were circulated within the wellbore and remains within the wellbore 
during testing and does not refer to the formation lluid that is located within the formation matrix 
or the formation fluid that is removed from the formation by action of the sampling tool or probe. 

After mud circulation stops, formation temperatures gradually build up. Assuming the 
heat flux rate is zero after mud circulation, by applying the superposition principle, the following 
equation can be obtained to express a one-dimensional radial temperature distribution during the 
period of temperature buildup (At). 

InHk , Inllk , 

By estimating the wellbore fluid temperature Tm and the original formation temperature 
Te, the value of Q can be calculated from equation (3). A calculated temperature T at a radial 
distance from the wellbore wall can be obtained versus lime from equation (4), an example of 
which is shown graphically in Figure 4, where curve 1 illustrates the temperature profile during 
mud circulation and curves 2-5 illustrate the temperature profiles at incremental times after 
circulation had stopped. 
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Fluid flow prediction 



The next step in the methodology of the present in\ ention invoK es ealculation of three- 
dimensional fluid How through spaee. 

Fluid withdrawal at the sink probe can be treated as a continuous point source with 
strength q, on a cylindrical internal boundary of radius r„ through which no temperature flux 
passes. For spherical fluid flow through space at a constant temperature caused by a continuous 
point source. Goode et al presented the pressure change at a point and time as follows: 



A/;(/) = 



a 



+ 



n, 



till 



(5) 



Where q, is an approximately constant withdrawal rate of a single phase liquid and (ro, 0„, 
Zo) are the cylindrical coordinates of the point source, and: 



a = r- + 1\; - 2/v-„ cos(6' - 0„ ) 
k. 



(PJUC, 



(PJLIC, 



(6) 
(7) 

(8) 



In order to obtain fluid flow equations, derivatives of pressure are taken with respect to r. 
0, and z. The pressure gradients along r, 8, and z may then be written as: 
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dp _ A /'-/■„ cos{ 0-6,, 
dr ~ B 2/7, 



B 

—J- exp 



V;r , Jb 

4jB 4t ) 



(9) 



dp _ A sin(6' - 6^,1 
de ~ B 21], 



T^cxpl 1+ - ^ ~e)ic 



' ^B' 



I ) 4^B 



(10) 



dp _ /[_^ 
d= ~ B Iq^ 



2V/ \ ( J 4Jb ■ 



f Jb ] 



(11) 



Where erfc(x) is the error function, and 



A = 



'I, 



(12) 



B 



a 



(13) 



Thus the velocity of a point in three-dimensional space at time t can be expressed as: 



// de r- 



(14) 



(15) 



V 



{') = - 



A_ dp 
~Jld^ 



(16) 



and 



A.V = vAf 



(17) 
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Where s is the distance of a space point moving in the direction of r, 0, z within the time 
period At and \' represents the velocity along the three directions. This is an analytic model, with 
which the location and velocity of a point in space can he calculated at any time. 

With this set of equations, the movement of fluid into the sink prohe at different times 
can be tracked, and hence the starting point can be located. Therefore, by dividing a small sphere 
around the point source into a family of streamlines, the starting point of the streamlines at a 
particular time can be located. The starting point of the fluid mo\ ement represents the location 
of that fluid at the start of the testing period. Figure 5 illustrates a portion of the solution of 
10 equations (5) through ( 17) at various drawdown times is shown in Figure 5. Once the locations 
of the starting points of the streamlines are known within the formation, the various locations can 
be referenced as a radial distance from the vvellbore wall. 

Due to the relatively short testing time involved with use of sampling tools, it is assumed 
that the temperalurc distribution within the formation is kept constant during testing. Based on 

15 the radial temperature distribution prediction as expressed in equation (4), the temperatures at 
these starting points can be calculated. The average value of the starting point temperatures 
represents a calculated temperature of the fluid at the starting point. The fluid flow equations are 
capable of estimating the time at which the fluid from each starting point will reach the probe. 
By knowing the temperature of the fluid at the starting point and the time at w hich this particular 

20 fluid will enter the probe, a calculated temperature of the fluid at the probe can be determined 
with respect to time. Although this is an analytic model, a computer program may be used to 
complete the tracking of fluid movement along the streamlines to the probe. 
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The solving of equations (9) through (17) provides a series of solutions of radial distance 
(location of the fluid) versus time, which are then combined with the temperature versus radius 
solution of equations (1) through (4) to provide a calculated temperature at the probe at a 
particular time. An example of a graphical representation of the calculated temperature versus 
5 time solution is shown as curve 2 in Figure 6. Varying the estimates of the wellbore fluid 
temperature Tm and the original formation temperature Te results in altering the calculated 
temperature versus time solution. 

Formation fluid measurements 

to A downhole probe is inserted into the wellbore where it can receive formation fluid 

samples and take temperature measurements of the formation fluid over the testing time period. 
Typically the downhole probe is engaged to the wellbore wall where it can receive formation 
fluid without contamination from the wellbore fluids, fhe formation fluid is withdrawn from the 
reservoir at a substantially uniform withdrawal rate. The formation fluid temperature 

t5 measurements are recorded and provide an observed formation fluid temperature at the wellbore 
versus time data. 



Identifying original formation temperature through optimization 



l o estimate the static formation temperature, the following objective function of the error 



20 E, between the observed and calculated temperatures is minimized: 




1 1 



(IX) 
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The wellbore fluid (mud) temperature Tm and initial formation temperature Te, are both 
adjustable parameters. The shape and slope of the calculated temperature plot is sensitive to the 
value of Tm, as well as the initial formation temperature Te. The observed temperature Tohs of 
the fluid entering the probe with respect to time can then be compared to the calculated 
temperature T,ai of the fluid entering the probe with respect to time. The error K, between these 
two temperature profiles is quantified in the solution of equation (18). yielding a measurement of 
the eiTor between them. Adjustments can be made to the estimates of wellbore fluid temperature 
Tm and original reservoir temperature Te to minimize this error. By systematically adjusting the 
wellbore fluid temperature Tm and original reservoir temperature Te values, the error between 
the calculated and measured temperatures can be minimized, resulting in an accurate prediction 
of the initial reservoir temperature Te. This process provides a systematic method of determining 
the static reservoir temperature, without requiring a prolonged temperature buildup period. 

Figure 7 is a schematic process diagram illustrating the methodology of one embodiment 
of the present invention. Box 1 contains the input data that is available, such as reservoir 
parameters, wellbore data and any pretest data that has been obtained. Box 2 represents the 
three-dimensional model used to analyze the formation surrounding the wellbore. The modeling 
of three-dimensional fluid flow through space allows a determination of formation fluid 
velocities at any point in the formation surrounding the wellbore, at any time. This technique is 
used to track the position of fluid at any time during the testing, and to determine the time at 
which fluid at a particular position will enter the wellbore, referred to as the time of flight. 
Streamlines of the fluid flow are generated and traced which enables the determination of a fluid 
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position in space as a starting point of fluid that enters the wellbore at a certain time. Box 3 
shows the parameter of time of flight location in space of formation fluid that w ill be used within 
further analysis. Bo.\ 4 contains the uipul data that is ax ailable relating to the thermal properties 
of the formation and wellbore. fhis data will include an estimate of the drilling fluid 
temperature within the wellbore and the static formation temperature. Bo.x 5 illustrates the 
temperature model that is used to convert the location in space data obtained from the flow 
model into a temperature field. The calculated temperatures for fluids having the same time of 
flight are averaged to give a calculated fluid temperature at various time of flights (i.e.. times at 
which particular formation fluids enter the wellbore), thereby giving an estimate of the fluid 
temperature entering the wellbore versus time. Box 6 shows the parameter of time of flight 
average temperature that gives a calculated fomiation fluid temperature entering the wellbore 
versus time that will be utilized for further analysis. Box 7 illustrates the comparison of the 
calculated temperature versus time data to the measured temperature versus time data that is 
obtained from the sample probe. The calculated formation fluid temperatures are based on the 
estimated temperature distribution prior to probe testing. The mud temperature in the wellbore at 
the end of wellbore fluid circulation and the initial formation temperature are both adjustable in 
the process of developing calculated temperature (entering the probe) versus time data. Box 8 
illustrates how the adjustable parameters of mud temperature and the initial reservoir temperature 
can be modified and input back into the temperature model shown as Box 5. Since the calculated 
temperature data is modified by adjusting the mud temperature and the initial reservoir 
temperature, by repetitive iterations, the error between the calculated and measured temperature 
data can be minimized, thereby providing an estimate of the static reser\ oir temperature. 
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Figure 4 graphically illustrates the temperature distribution after wellbore fluid circulation 
has ceased. Temperature distributions along a radial distance from the wellbore at several build 
up stages are shown. After wellbore fluid circulation stops, the temperature funnel rapidly builds 
up toward a uniform temperature field, but the buildup becomes increasingly slower, which 
5 indicates heat flux decrease with time. 

Figure 5 graphically illustrates the calculated three-dimensional fluid flow profiles at 
various drawdown times. Due to permeability anisotropy, the flow profiles are elliptic and 
flatten with time. 

Figure 6 graphically illustrates the calculated and measured temperatures at the sink probe 
10 versus elapsed time after wellbore fluid circulation has ceased. The relationship of the 
temperature measured at the MDT probe and the calculated temperature versus time is shown for 
the application example. Good matching results of the calculated temperature plot with the 
measured temperature plot were obtained for most of the withdrawal time. The early data where 
the measured data is significantly higher than the calculated data indicates that the initial MDT 
15 flow line temperature was higher than the fluid flowing into the tool at that time, however, the 
measured temperature rapidly decreased to the actual flowing fluid temperature. Prior to the 
application example testing, the MDT tool had tested deeper intervals (having higher 
temperatures) during the same run prior to testing the point of interest, which may explain the 
deviation in the early data. The good match with field data indicates that this approach is 
20 feasible and reliable. 



Page 19 of 28 



PATENT 
20.2787 

One particular embodiment of the present invention involves the estimating of the static 
formation temperature and the mud temperature within the wellbore after mud circulation has 
stopped. Using these estimates along with physical properties of the reservoir, circulation times 
and buildup times (length of time after circulation has stopped), the one-dimensional radial heat 
5 flux equations of equations ( 1 ) through (4) can be solved. These solutions create a formation 
temperature profile versus the radial distance from the wellbore at various buildup times, such as 
shown graphically in Figure 4. Using an estimated formation fluid withdrawal rate at the 
wellbore by a sample probe, a three-dimensional fluid flow model of the reservoir can be 
calculated from equations (5) through (17), establishing a formation fluid location profile at 

10 various drawdown (sampling/testing) times, such as shown graphically in Figure 5. Knowing the 
drawdown time and the buildup time, the starting points of the fluid movements can be located. 
These starting point locations can be related to the temperature profile versus radial distance plot 
in Figure 4, determining the temperatures for each of these starting locations. These 
temperatures can be averaged to predict the temperature of the fluid sample that will comprise 

15 the combination of fluid from these locations. 

By predicting the time that this particular fluid will enter the sample probe from the fluid 
flow model, and the predicted temperature that this particular fluid will have as it enters the 
sample probe, a predicted temperature versus time profile for the formation fluid that is 
withdrawn from the formation by the sample probe can be developed. By inserting a sample 
20 probe within the wellbore and engaging it with the wellbore wall, formation fluid can be 
withdrawn from the formation. The fluid is ty{")ically withdrawn at a substantially known rate 
and at a substantially steady rate to reduce variations in fluid flows through the probe that can 
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alter the test data. If the fluid withdrawal rate is substantially different than the estimated fluid 
withdrawal rate used in the fluid How model, the fluid How may need to be recaleulated with the 
actual tluid flow rate. The temperature of the formation fluid passing through the probe at the 
wellbore is measured, which creates an obser\ ed formation fluid temperature versus time profile. 

The error between the observed formation fluid temperature at the wellbore versus time 
profile and the calculated formation fluid temperature at the wellbore versus time profile can be 
quantified mathematically, such as by equation (18). Since the calculated formation fluid 
temperature at the wellbore \ ersus time profile is dependant on the estimates of the static 
formation temperature and the mud temperature within the wellbore, these estimates can be 
10 modified and a revised calculated formation fluid temperature at the wellbore versus time profile 
generated. 

The predicted static formation temperature can be obtained by repeating the iteration of 
modifying the estimated static formation fluid temperature and wellbore fluid temperature 
variables, and generating a revised calculated fomiation fluid temperature at the wellbore versus 
\5 time profile, until the error between the calculated formation fluid temperature at the wellbore 
versus time profile and the observed formation fluid temperature at the wellbore versus time 
profile is minimized. 
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